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Abstract. Muons with a high transverse momentum 
(Pt) are produced in cosmic ray air showers via 
semileptonic decay of heavy quarks and the decay 
of high p T kaons and pions. These high p T muons 
have a large lateral separation from the shower core 
muon bundle. IceCube is well suited for the detection 
of high p T muons. The surface shower array can 
determine the energy, core location and direction 
of the cosmic ray air shower while the in-ice array 
can reconstruct the energy and direction of the high 
p T muon. This makes it possible to measure the 
decoherence function (lateral separation spectrum) at 
distances greater than 150 meters. The muon p T can 
be determined from the muon energy (measured by 
dE/dx) and the lateral separation. The high p T muon 
spectrum may also be calculated in a perturbative 
QCD framework; this spectrum is sensitive to the 
cosmic-ray composition. 
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I. Introduction 

The composition of cosmic rays with energies above 
10 6 GeV is not well known. At these energies, the 
flux of cosmic rays is so low that direct detection 
with balloon or satellite-borne experiments is no longer 
feasible and indirect measurement with larger ground 
arrays must be used. These arrays measure the electronic 
and hadronic components of a cosmic ray air shower, 
and must rely on phenomemological interaction models 
to relate observables like muon and electron density 
to primary composition. These interaction models are 
based on measurements made at accelerators that reach 
a maximum energy roughly equivalent to a 10 6 GeV 
proton cosmic ray (TJ. Extrapolation to the energy of 
the detected cosmic rays leads to uncertainties in the 
composition of cosmic rays at high energies. An alter- 
nate method of determining the composition is to use 
muons with a high transverse momentum (px) 0. 

At transverse momenta on the order of a few GeV/c, 
the muon px spectrum can be calculated using pertur- 
bative QCD (pQCD). Such calculations have been made 
for RHIC and the Tevatron, and the data is in quite good 
agreement with modern fixed order plus next-to-leading 
log calculations [3]. These experimental studies give us 
some confidence in pQCD calculations for air showers. 

Most of the high-energy muons that are visible in 
IceCube are produced from collisions where a high- 



energy (large Bjorken— x) parton interacts with a low 
Bjorken— x parton in a nitrogen or oxygen target. These 
collisions will produce heavy (charmed and bottom) 
quarks and also jets from high px patterns; the jets will 
fragment into pions and kaons. Pions and kaons produce 
"conventional" muons that have a soft spectrum roughly 
proportional to E -3 ' [4] and typicallly have a low px- In 
contrast, charm and bottom quarks are produced early in 
the shower. The resulting muons ("prompt" muons) have 
a harder spectrum, a higher px, and are the dominant 
source of atmospheric muons above 10 5 GeV Q. 

If these particles are produced in the forward direction 
(where they will be seen by ground based detectors), 
they can carry a significant fraction of the energy of 
the incident nucleon. The muon energy and px can be 
related to the energy of the partons that make up the 
incident cosmic ray. For example, a 10 16 eV proton 
has a maximum parton energy of 10 16 eV, while the 
maximum parton energy for a 10 16 eV iron nucleus 
(with A=56) is 1.8 x 10 14 eV. These two cases have 
very different kinematic limits for high px, high-energy 
muon production, so measurements of high px muon 
spectra are sensitive to the cosmic-ray composition. 

High px muons constitute a small fraction of the 
muons visible in IceCube; the typical muon px is of 
order a few hundred MeV/c. In this domain of "soft 
physics," the coupling constants are very large, and 
pQCD calculations are no longer reliable. One is forced 
to rely on phenomenological models; different hadronic 
interaction models give rather different results on com- 
position o. 

IceCube, a kilometer-scale neutrino telescope, is well 
suited for the detection of muons with px& above a few 
GeV/c [2|. When completed in 2011, it will consist of a 
1 km 3 array of optical sensors (digital optical modules, 
or DOMs) buried deep in the ice of the South Pole and a 
1 km 2 surface air shower array called IceTop. IceTop has 
an energy threshold of 300 TeV and can reconstruct the 
direction of showers with energies above 1 PeV within 
~1.5° and locate the shower core with an accuracy 
of ~10 m [7|. The in-ice DOMs (here referred to as 
IceCube) are buried in the ice 1450 m under IceTop in 
kilometer-long strings of 60 DOMs with an intra-DOM 
spacing of 17 m. IceCube can reconstruct high quality 
tracks of high energy muons with sub-degree accuracy. 
IceTop can measure the energy, position, and direction of 
the shower and IceCube can do the same for the muons. 
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These values can be used to calculate the pt'- 



where is the energy of the high pt muon, d is its 
lateral separation and h is the interaction height of the 
shower, here taken as an average value of 25 km. The 
interaction height loosely depends on the composition 
and a full treatment of this is planned in the future. Tak- 
ing 150 m, 25 meters more than the separation between 
strings of DOMs in IceCube, as a rough threshold for 
the two-track resolution distance of the high pt muon 
from the shower core gives a minimum resolvable px of 
6 GeV/c for a 1 TeV muon. 

Multiple scattering and magnetic fields can deflect the 
muons as they travel to the IceCube detector, but this 
is only equivalent to a few hundred MeV worth of pr 
and is not a strong effect in this analysis. The high pt 
muon events are near-vertical which gives them a short 
slant depth. In order for muons to reach the IceCube 
detector they must have an energy of about 500 GeV 
at the surface of the Earth. These higher energy muons 
are deflected less by multiple scattering and the Earth's 
magnetic field. 

The combined acceptance for cosmic ray air show- 
ers which pass through both IceTop and IceCube is 
0.3 km 2 sr for the full 80-string IceCube array (8). By 
the end of the austral summer of 2006/2007, 22 IceCube 
strings and 26 IceTop tanks had been deployed. The 
combined acceptance for showers that pass through both 
IceTop-26 and IceCube-22 is 0.09 km 2 sr. While this ac- 
ceptance is too small to expect enough events to generate 
a px spectrum, it offers an excellent opportunity to test 
reconstruction and background rejection techniques. 

II. Previous Measurement of High p T Muons 

MACRO, an underground muon detector, has previ- 
ously measured the separation between muons in air 
showers with energies ranging roughly from 10 4 GeV to 
10 6 GeV [9|. MACRO measured muon separations out 
to a distance of about 65 meters. The average separation 
between muons was on the order of 10 m, with 90% of 
the muons found with a separation of less than 20 m [9]. 
MACRO simulated air showers and studied the muon 
pair separations as a function of the px of their parent 
mesons. They verified the linear relationship between px 
and separation shown in Eq. 03 (with a small offset due 
to multiple scattering of the muons) out to momenta up 
to 1.2 GeV/c. 

III. High p T Rate Estimations 

The decay of charm will produce roughly 10 6 prompt 
muons/year with energies in excess of 1 TeV inside the 
0.3 km 2 sr combined acceptance of the full IceCube array 
flOl . Different calculations agree well in overall rate in 
IceCube, but, at high (1 PeV) energies, different parton 
distribution functions lead to rates that differ by a factor 
of 3 HQ. 



TABLE I 

Estimated number of events from charm above different 
p T thresholds in 1 year with the 80-string IceCube array 



Pt [GeV/c] 


Separation [m] 


Number 


>6 


150 


~500 


>8 


200 


~200 


>16 


400 


~5 



We can roughly estimate the fraction of muons that 
are produced at high pt using PYTHIA pp simulations 
conducted for the ALICE muon spectrometer at a center 
of mass energy of 14 TeV, which indicate that roughly 
0.6% of these events will have a pt of at least 6 GeV/c 
lfl2l . These simulations are done at a higher center of 
mass energy than the bulk of the IceCube events, and 
also were for muons produced at mid-rapidity (and, 
therefore, higher projectile x values and lower target x 
values than in IceCube). However, they should be ade- 
quate for rough estimates. These estimated rates of high 
Pt muons are further reduced by approximately a factor 
of 10 by requiring that the high pt muon be produced in 
coincidence with a shower that triggers IceTop, leaving 
a rough expectation of ~500 events per year with a 
Pt greater than 6 GeV/c in the 80-string configuration 
of IceCube. The estimated number of events above a 
given pt from the decay of charm is given in Table [I] 
neglecting the uncertainties mentioned above. At higher 
Pt, bottom production may also contribute significantly 
to muon production. 

The rate of 1 TeV muons from conventional flux is 
expected to be about 10 9 events/year in the full IceCube 
configuration. The vast majority of these muons will 
have a low pt- Based on the pt spectrum measured 
from pions produced in 200 GeV center of mass energy 
pp collisions by the PHENIX collaboration, the number 
of events expected with pt > 6 GeV/c is estimated to 
be 1 in 6 x 10 6 lfl3ll . Requiring a pt > 6 GeV/c and 
an IceTop trigger leaves roughly 20 events/year in the 
full IceCube configuration from the conventional flux of 
muons. 

IV. Reconstruction Methods 

High pt muons will appear as a separate track coinci- 
dent in time and parallel with the track from the central 
core of low pt muons. Generally the bundle of low pt 
muons leaves more light in the detector than the high 
Pt muon. 

Current reconstruction algorithms in IceCube are de- 
signed to reconstruct single tracks. In order to recon- 
struct these double-track events the activated DOMs (i.e. 
DOMs that detect at least one photon) are split into 
groups using a k-means clustering algorithm [TBI . The 
first group is the muon bundle and the second group 
is the high pt muon. Each group is then reconstructed 
with a maximum likelihood method that takes into 
account the scattering of light in ice. After this initial 
reconstruction, the groups are resplit according to their 
time residual relative to the muon bundle reconstruction 
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Fig. 1. Zenith angle resolution of the high px reconstruction algo- 
rithm. The sigma are the results of Gaussian fits to the distributions. 



and re-fit. The first splitting forces the activated DOMs 
into two groups, which is not correct for the background 
events which don't have a high pr muon and generate 
only a single shower in the array. Splitting the activated 
DOMs a second time according to their time residual 
allows for the possibility for all the activated DOMs to 
end up in a single group. FigureQ]shows the performance 
of this clustering algorithm. The zenith angle resolution 
for groups determined using the true simulation infor- 
mation (black, solid lines) is compared to the resolution 
for groups determined using the clustering algorithm 
(red, dashed lines) for the muon bundle (top) and high 
Pt muon (bottom). Roughly 20% of the events fail to 
reconstruct because there are not enough DOMs in one 
of the groups. 

These reconstruction algorithms achieve a zenith an- 
gle resolution of 2.7° for the muon bundle and 4.1° for 
a high pt muon separated by 400 m. The resolution is 
worse for the high pr muon because fewer DOMs are 
activated. While high pt muons with a greater separation 
are much easier to resolve with the two track algorithm, 
they also tend to be lower energy (see Eq. [TJ and activate 
fewer DOMs. The average number of DOMs activated 
by the high pt muon is 12, compared to 50 for the muon 
bundle. Additionally, because high pr muons have less 
activated DOMs, a DOM activated by the muon bundle 
that is incorrectly placed in the high pr group has a 
much larger effect on the reconstruction of the high pt 
track (it can also degrade the bundle resolution, but to a 
lesser extent). These factors lead to a poorer resolution 
in the high pt track direction. 

V. Signal and Background Separation 

The interaction models in the air shower simulation 
program CORSIKA [15| do not include many of the 
processes that can generate a high pt muon, with the 
exception of the modified DPMJET model discussed in 
fl6l . Also, since high pt muons will occur in only a 
fraction of simulated showers making simulation very 
time intensive, a toy model based on CORSIKA proton 



showers was used to model the signal. A single muon is 
inserted into an existing CORSIKA event containing a 
muon bundle from an air shower. This modified shower 
is then run through the standard IceCube propagation, 
detector simulation, and reconstruction routines. Simu- 
lations insert a muon with energy of 1 TeV separated 
100, 150, 200, and 400 m from the shower core. This 
gave a clean sample ideal for the development of re- 
construction routines optimized to identify simultaneous 
parallel tracks inside the detector. 

Cosmic rays air showers that do not generate a high 
Pt muon are considered a background to this search. 
Since they generate only a single shower in the array, 
these events are mostly eliminated by requiring there be 
two well-reconstructed tracks in the IceCube detector. 
These single showers are well-reconstructed by a single 
track hypothesis, while the high pt muon events are not. 
Figure|2]shows the negative log of the reduced likelihood 
of a single track reconstruction for single showers, and 
showers with an inserted 4, 8, and 16 GeV/c pr, 1 TeV 
muon (separation of 100 m, 200 m, and 400 m from 
the shower core, assuming an average interaction height 
of 25 km). Well-reconstructed events have a lower value 
on this plot. For large separations, this variable separates 
single showers from showers which contain a high pt 
muon. When the separation between the high pt muon 
and the shower core drops below the interstring distance 
(the blue, dot-dashed line in Figure [2j, it is no longer 
possible to cleanly resolve the high pt muon from the 
shower core and the event looks very similar to a single 
shower. 



= 0.25 



£ 0.2 
< 



4 GeV/c Muon 
8 GeV/c Muon 
16 GeV/c Muon 
Single CR Shower 




7 7.5 8 8.5 9 9.5 10 
-Log (Single Track Reduced Likelihood) 



Fig. 2. Negative log of the reduced likelihood for the single track 
reconstruction 



The IceCube 22-string configuration is large enough 
that the rate of simultaneous events from cosmic rays 
is significant. Muon bundles from two (or more) air 
showers can strike the array within the 10 fjis event 
window, producing two separated tracks. These so-called 
double-coincident events are the dominant background 
for air showers with high pt muons. Since these double- 
coincident events are uncorrelated in direction and time, 
requiring that both reconstructed tracks be parallel and 
occur within 1 (is can eliminate most of these events. 
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Due to misreconstructions some events will survive 
this selection, and criteria to eliminate these misre- 
constructed events are being developed. However, an 
irreducible background remains from double-coincident 
events that happen to come from roughly the same 
direction and time. A conservative estimate of the num- 
ber of these events is comparable to the number of 
expected high pt events. A full simulation is under way 
to determine the rate of double-coincident cosmic ray 
events and to develop methods to distinguish high px 
events from this background. 
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Fig. 3. Candidate shower with a high px muon. The colored circles 
show DOMs that are activated by light. The color and size of the 
circle corresponds to the time (red being earliest) and magnitude of 
the signal, respectively and the white dots indicate DOMs that are 
not activated in the event. The red lines show the reconstructed tracks 
(with red circles denoting the arbitrary vertex of the reconstruction). 

After applying selection criteria to reduce the cosmic 
ray background, a 10% sample of the IceCube data was 
scanned by eye for high pt muon candidates. Several 
events were found and Figure [3] shows a representative 
event. The two tracks occur within 400 ns of each other 
and the space angle between the two reconstructions is 
5.6°. 

VI. Horizontal Events 

One possible method to avoid the backgrounds from 
single and double-coincident cosmic ray air showers is to 
search for events that come from the horizontal direction. 
At a zenith angle of 60° the slant depth through the 
ice is about 3 km. Only muons with energies in excess 
of several TeV will reach IceCube ifTTl and the event 
multiplicity is greatly reduced. Because the intra-DOM 
spacing is 17 m, compared to the intra-string spacing 
of 125 m relevant for vertical events, the two-track 
resolution should improve. A disadvantage to searching 
in the horizontal direction is the increase in the deflection 
from the magnetic field and multiple scattering with the 
slant depth, which could lead to a greater spread in the 
muon bundle and create event topologies which mimic 
a high pt muon event. Nevertheless, the advantage of 
background suppression is strong. In addition to high 
Pt muons in cosmic rays, there are a number of other 
processes (such as the decay of the supersymmetric stau 
010) that can produce horizontal parallel tracks, making 
the horizon an interesting direction. 



Fig. 4. Candidate double track near-horizontal event 

A search was conducted in 10% of the IceCube data 
from 2007 for these horizontal events using simple topo- 
logical selection criteria. Events which reconstructed 
within 30° of horizontal were selected and searched for 
a double track topology. Several interesting candidate 
events were found, one is shown in Figure [4] 

VII. Conclusions 

IceCube is large enough that the study of high pt 
muons has become viable. Resolution of tracks separated 
by at least 150 meters from the shower core will allow 
identification of muons with px& of at least 6 GeV/c. A 
two-track reconstruction algorithm has been developed 
and selection criteria for identification of air showers 
with high pt muons is under development. A search for 
horizontal double tracks is also under way. The rate of 
high pt muon production is sensitive to the composition 
of the cosmic rays and offers an alternative to existing 
composition studies. 

This work is made possible by support from the NSF 
and the DOE. 
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